Irreversible nucleation in molecular beam epitaxy: Prom theory to experiments 
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Recently, the nucleation rate on top of a terrace during the irreversible growth of a crystal surface 
by MBE has been determined exactly. In this paper we go beyond the standard model usually em- 
ployed to study the nucleation process, and we analyze the qualitative and quantitative consequences 
of two important additional physical ingredients: the nonuniformity of the Ehrlich-Schwoebel barrier 
at the step-edge, because of the existence of kinks, and the steering effects, due to the interaction 
between the atoms of the flux and the substrate. We apply our results to typical experiments of 
second layer nucleation. 
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I. INTRODUCTION 

Nucleation is a key process during the growth of a crys- 
tal by Molecular Beam Epitaxy (MBE), when the sub- 
strate is oriented along a high symmetry direction. In 
that case, freshly deposited atoms diffuse on the surface 
until they meet other adatoms (nucleation process) or the 
step of a growing terrace (aggregation process). Layers 
are completed through island coalescence and the filling 
of vacancies. 

Nucleation is irreversible if a dimer (i.e., a nucleus 
formed by the meeting of two adatoms) is thermally sta- 
ble. This condition depends on the two most important 
external parameters: the temperature T and the inten- 
sity F of the flux. It holds if the dissociation time of 
a dimer (which grows with decreasing T) is larger than 
the time (decreasing with F) required to aggregate more 
adatoms to the dimer and to stabilize it. 

As a matter of fact, nucleation appears to be irre- 
versible, in a range of temperatures experimentally ai 
cessible, for several experiniental systems as PtyPt(lll),l 
Ag/Ag(100),ElFe/Fe(100),y and Ag/Ag(lll).y 

Except for the submonolayer regime, where adatoms 
diffuse on the substrate, nucleations occur an terraces, 
and most of them take place on 'top terraces',El which are 
defined as those surrounded by a single closed descending 

In some recent paper J^'i we have studied the nucleation 
process on top of a terrace: we have evaluated the number 
Lv of nucleation events per unit time (nucleation rate) and 
their spatial distribution. The total nucleation rate u is 
of great importance in several respects: it determines 
the typical distance between nucleation centers in the 
submonolayer regime, the so-called diffusion length £o 
(see Ref. |^); it allows to extract the additional energy 
barrier for interlayer diffusion [Ehrlich-Schwoebel (E~ 
barrier] from an experiment of second layer nucleation;! 
it is the main factor determining the stable or unstable 
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character of growth during its first stagesB 

In Sec. H of this paper we report the formula giving 
the exact Yf-H^^ ^f the nucleation rate uj that was derived 
previouslyE^S Such a formula depends on the quantity 
W, the probability that two adatoms deposited simulta- 
neously actually meet. Here we evaluate W explicitly for 
all values of the terrace size L and of the so-called ES 
length, — 'p' ~ which measures the difference be- 
tween the intra-layer hopping rate v and the inter-layer 
hopping rate v' . The resulting formula is extremely ac- 
curate. In this way uj is written in full generality in terms 
only of L, F, v, v' and some geometrical constants that 
are given explicitly for the most relevant cases. We also 
provide the approximate formulas which are correct in 
some limiting regimes. 

Later we discuss, both from a qualitative and a quanti- 
tative point of view, two issues which are believed to be 
important experimentally and that are not included in 
the 'standard' models used for studying nucleation: the 
effects of nonuniform interlayer barriers and of steering 
phenomena in the deposition flux. 

The 'standard' model for nucleation assumes the pres- 
ence of an interlayer barrier which is uniform along the 
step surrounding the top terrace. However, even for com- 
pact terraces, kinks are unavoidably present along, the 
step and they are preferential sites for descendingja be- 
cause the ES barrier is expected to be reduced by the 
increased number of neighbours. The 'effective' barrier 
ff is therefore a combination of the barrier Ai?fc at 
kinks and the higher barrier Ai?o, felt by an adatom 
along straight steps. Even more importantly, since the 
number of kinks is temperature and coverage dependent, 
Ai?„tt depends on T and coverage, even if AiJ^ and Ai?o 
do not. This issue is treated in Sec. [II. 

The 'standard' model also assumes a uniform flux 
F impingipg on the terrace. However, it has been 
shownliErEJilj that the attractive interaction between the 
incoming atom and the atoms incorporated in the grow- 
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ing surface increases the number of particles landing on a 
top terrace. This may clearly affect the nucleation rate. 
The relevance of this effect strongly depends on the ki- 
netic energy of the incoming particles, the strength of 
the interaction and the angle of deposition. A minimal 
model, where all these factors determin e a s ingle adimen- 
sional parameter e, is discussed in Sec. pV} 

In Sec. we apply our results to typical experiments of 
second layer nucleation: an exppjiment on Pt/Pt(lll) at 
different CO partial pressures,tj and an experiment on 
Ag/Ag(lll) at different temperaturesJj In this way we 
are able to discuss the relevance of nonuniform barriers 
and steering effects in two real cases. 



Concluding remarks are presented in Sec. Vl 



II. THE NUCLEATION RATE 

In Refs. ^ and ^ we have evaluated exactly the nucle- 
ation rate w for a terrace of linear size L. This quantity is 
adimensional and represents the number of atoms along 
the edge, for a polygon-shaped island, or along the ra- 
dius, for a circular island. The total number of atoms in 
the island \s A = jL'^- The most general expression for 
the nucleation rate is: 



u^T ^-^ 



(1) 



where J- is the number of atoms arriving on the terrace 
per unit time, r^ep and r^-es a-re the deposition and res- 
idence time, respectively, and W is the probability that 
two adatoms, deposited simultaneously, meet before leav- 
ing the terrace. 

For a uniform flux, !F = FA, and its inverse is the 
average time between deposition events, the deposition 
time, Tdep — The residence time Tres is the average 

time spent by a single particle on the terrace before de- 
scending from it. It depends on the size L of the terrace 
and on the strength of the additional ES barrier, which 
can be quantified via the (adimensional) ES length, ^es = 
■p- — 1. The rates v, v' for intra and inter-layer hopping — 
whose precise definition is given in App. ^ — are gener- 
ally written as a prefactor times an Arrhenius factor, 
u = voeM-Ed/KT), v' = j/^exp[-(i;d + l\E)lk^T\. 
In the simple case whece vq and fg are supposed the 
same, the ES length iiil 4s = <sx^{l^E j k^T) In 
terms of L and £es , the residence time is given by 



(a4s + m^l^ : 



(2) 



which allows to distinguish in an easy way the 'weak bar- 
rier' regime (4s ^ Tres — PL'^/v) from the 'strong 
barrier' one (4s ^ r^es — oiLi^s/i^ = aL/v'). The 
numerical factors a and as well as 7, depend on the 
shape of the terrace and on the symmetry of the under- 
lying lattice. Their values for some relevant cases have 
been determined numerically or analytically and are re- 
ported in Table |[ In App. ^ we determine analytically 
the geometrical parameters relevant for a circular terrace. 



It is worth noting that the parameter a can easily be 
written in terms of 7 and other geometrical factors. As 
a matter of fact, in the limit of large ES barriers, the 
residence time has two equivalent expressions. The first 
comes from Eq. (||), Tres = aL/v' . The second is found 
by considering that the escape rate from the terrace, r^^, 
is given by the probability of finding an adatom on an 
edge site, times the fraction of hops leading from an edge 
site to the lower terrace, times the interlayer hopping 
rate. The fraction of jumps leading to the lower terrace is 
equal to Az/z, where z is the coordination number of the 
lattice, and Az is the number of missing neighbours for 
an edge site. In the strong barrier regime the probability 
distribution of adatoms on the terrace is uniform and 
hence the probability of finding an adatom on an edge 
site is V I A, where V = rjL is the perimeter of the terrace 
in units of lattice sites. Summing up, we can write 



A z 
V~Kzv 



1 



7z 



TyAz v' 



[4s » L] (3) 



so that a — / (ril^z). The quantities rj and Az are 
reported in Table | (see also App. 

The quantity W in Eq. (|l]), the probability that two 
atoms both on the terrace actually meet, depends on the 
initial spatial distributions for the two atoms, but this d|&=. 
pendence is very weak and can in practice be neglected.E£l 
The probability W does not change much even if one 
atom is taken as immobilcEl In such a case, W can easily 
be written in terms of properties of a single random walk 



W 



A 



(4) 



where N^^^^L, l^^^) is the number of distinct sites visited 
by a single particle before leaving the terrace and clearly 
depends on £^3 and L. For extremely strong barriers, the 
atom is able to visit all sites, N^^^ A, so that W ^ rL. 
In the opposite limit of very weak barriers, it is knownliZI 
that N^UL,0) = NoL^/\n{L/Lo) (so that W « 1/lni). 
An approximate formula interpolating between the two 
limits is given in Ref. ^. Here we improve on that esti- 
mate by providing an analytic expression that reproduces 
with very good accuracy the value of W obtained by nu- 
merical simulations for all values of £es- The formula 
is 
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W 



NoL\n{L/Li) 
HL/Lo) 



Nii^s + lLln{L/Li) 



(5) 



The values of Nq, Lq, Ni, Li depend on the shape of the 
terrace and on the symmetry of the underlying lattice, 
as a, (3 and 7 do. Their values are given for relevant ter- 
race shapes and lattice types in Table |. The derivation 
of Eq. (||) is given in App. ^ and its accuracy is demon- 
strated in Fig. |l]. Notice that no parameter is fitted. 

By inserting Eq. (||) and Eq. (|) into Eq. (|l|) we ob- 
tain the general fully explicit formula for the nucleation 
rate. Eq. (|l|) is based on the sole hypothesis that the 
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TABLE I: Numerical values of the parameters q and /3 [ap- 
pearing in Tre.s, Eq. (^)], 7 (appearing in 77 (appearing in 
V), Az [appearing in Tr^s, Eq. (|)], No,Lo,Nr,L-i_ [appearing 
in W, Eq. dsl)]. They are given for different shapes of the 
terrace and for two types of lattice: the square and the tri- 
angular one. They are typical of the (100) and (111) face of 
a cubic lattice, respectively. For triangular (A), square (□), 
and hexagonal (o) islands, L is the edge of the polygon. For 
circular (Q) terrace, L is the radius. The values of a are 
computed using the formula a = ^z/{r]A.z) [see Eq. (^] and 
agree with numerical results. The values of A^o, Lq, N\ and Li 
are determined by fitting numerical results for Ndis{L) and 
Af^f (L) to their analytical expressions. 



by the expression: 



C\z No Lo Ai Li 



(100) □ 1 0.14 1 4 1 0.2 1.0 1.8 0.5 

O vr/2 0.5 TT 4^2 0.7 0.5 3.1 0.3 

(111) A 1/2 0.05 1/2 3 2 0.08 2.3 0.9 1.6 

O 3/2 0.4 3 6 2 0.7 0.5 3.0 0.4 

O vr/2 0.5 4^3 2 0.7 1.0 3.8 0.2 




FIG. 1: Comparison of the exact numerical values of W 
(points, see Ref. ^), with the approximate analytic formula 
(line) given in Eq. (p|), for a square lattice and a square ter- 
race of size L = 20. 



deposition time Tdep is much larger than the so-called 
'traversal time' Ttr = j3L'^ /v, defined as the average time 
an atom needs to reach the edge of the terrace (note 
that Ttr — Tresi in the reginae £es 

< L)B The condi- 
tion Tdep 3> Ttr is cquivalentu to the condition 3> 1. 
A quantitative evaluation for the experimental systems 
discussed in Sec. gives in all cases Ttr /Tdep < lO^'"^. 

In most cases, also the residence time is much smaller 
than the deposition time. For example, for the experi- 
mental systems studied in Sec. ^we have r^es/Tdeplpt < 
5 X lO^'^ and Tres/Tdep\Ag < 10~^. If the condition 
Tres ^ Tdep is Satisfied, the nucleation rate is just given 



' res 
I2 



W [Tres < Tdep] 



f2 



(6) 



(7) 



where W is given in Eq. (^). In the limit of strong bar- 
riers, if we papproximate W by one, we have the simpler 
expression:E3 



ui = q;7 



[Tres < Tdep and — > i] (8) 



III. THE EFFECT OF NON UNIFORM 
BARRIERS 

An adatom leaves the terrace because it reaches an 
edge site and it jumps down with a hopping rate v' : in the 
usual model for nucleation, v' is supposed to be the same 
along the step edge. However, its value may actually be 
non uniform: a growing step always has some degree of 
roughness, which makes the descent preferable from cer- 
tain edge sites, and toward certain directions. For an. 
fcc(lll) surface, the picture is even more complicated,^ 
because steps can be of two different types (A-step and 
B-step) and they can both be present around the same 
terrace. In this case their barriers are different even if 
both types of steps are straight. Therefore, if the stan- 
dard model is adopted, the interlayer hopping rate should 
be understood as an effective one, valid at well defined 
temperature and growth conditions. Any factor able to 
affect the step morphology may affect v' as well. The 'ef- 
fective' character of v' makes its separation in a prefactor 
and an exponential factor arbitrary. For this reason we 
arc going to assume Vq = vq. 

The case of a generic distribution of ES barriers can 
not be treated analytically, because it is not possible 
to find the general expressions for Tres and W . In the 
following we are assuming to be in the 'strong barrier' 
regime, where W = 1 and an analytical derivation of Tres 
is indeed possible. In this regime, the adatom density is 
uniform and Eq. (^ is easily generalized to non-uniform 
barriers: 



A 



1 



V (Az) {v') 



(9) 



Here (Az) is the number of directions leading to a hop 



and 



is the 



downwards (averaged over edge sites), 
average hopping rate. 

Let us consider a simplified 'kink- model': the ES bar- 
rier is everywhere equal to the high value Ai?o, except 
for special sites (and paths), where the additional bar- 
rier takes the smaller value AEk- In other words, the 
distribution is bimodal 



AEk with probability Ck 
AEq with probability (1 — Ck). 



(10) 



4 



Eq. is valid as long as even the smallest of the bar- 
riers is sufficiently large to ensure a uniform probability 
of finding the adatom on the terrace. Hence the criterion 
for its validity is that the smallest ES length is much 
larger than the linear size of the terrace, 



ES ifc 



(11) 



Provided this is true, the effective barrier is determined 
by the relation 



(12) 



whose general solution is 



A£;,., = ASfc + fceTln 



Cfe 



(13) 



-fc„nn 



Cfe 



1 exp[-{AEo- AEk)/kBT] 



In the case when the two barriers are practically the 
same, Eq. (|l^ obviously gives AE^a « AEk w AEq. 
In the more interesting opposite limit, when AEq — AEk 
is larger than kBT\n{l/ck — l), the last term on the right- 
hand-side is negligible, and the effective barrier has the 
simplified expression, 



AE,!f ^ AEk + fcerin 



Cfe 



(14) 



Eq. (|Ij) has a transparent physical meaning: adatoms 
leave the terrace at kink sites only, but they feel an effec- 
tive barrier larger than AEk, because of the finite con- 
centration of kinks. If T = 300 K and Cfe ~ 0.1, the 
barrier increase due to this effect is 0.06 cV. 



What happens when the condition (11) does not hold 
because AEk is not large enough? In that case the prob- 
ability to find the adatom at a kink site is suppressed and 
the use of Eqs. (H) and (14) underestimates r^es- This is 
clearly shown in Fig. ^ where we compare, for a square 
terrace on a square lattice, numerical results (symbols) 
for Tres (in units l/v = 1) with the analytical approx- 
imation (lines) derived from Eqs. (||) and (|l^. In the 
case AEq = oo and AEk — the analytical formula 
gives a residence time (dashed line) that is well below 
the numerical results (diamonds). In the opposite case 
of AEo = oo and (4s)fe = exp(A£;fe/A:Br) - 1 = 99, the 
condition (O) is fulfilled: numerical (circles) and analyt- 
ical results (solid line) agree. 

Eq. (|l^) clearly applies not only to the case of kinks, 
but holds whenever there are two energy barriers for in- 
terlayer transport, one of which is significantly larger: 
For example for an hexagonal terrace on a (111) surface, 
surrounded by three A-steps and three B-steps. If AEa,b 
are the ES barriers at the two different steps, the effective 
barriet.is equal to the smaller one plus a small correction 
termJij 



AE^a = miTi{AEA, AEb} + k^T In 2 



(15) 




FIG. 2: The residence time (in units Xjv = 1) for an adatom 
on a square terrace of size L = 40, as a function of the in- 
verse of the kink concentration (log- log scale). Ai5o = oo 
in all cases, while A_Efc corresponds to (^Es)fc = for dia- 
monds and to (£Es)fc = 99 for circles. Symbols refer to the 
exact numerical calculation of Tres and lines to the analytical 
approximation, Eqs. (PlH). 



IV. STEERING EFFECTS 

The flux of incoming atoms is usually supposed to be 
uniform. This hypothesis is correct down to a distance 
of a few nanometers from the surface, but it breaks down 
at smaller jiistances because of the adatom-substrate 
interaction.E£l This interaction is attractive, so tha|tj|Ljn- 
creases the effective flux landing on a top terrace. Ot3 

The effect is illustrated in Fig. || for normal deposition 
on a circular terrace. Because of the attractive interac- 
tion with the terrace, the deposited atom may deviate 
from its rectilinear trajectory t\ and land on the terrace. 
The incoming flux JF, equal to FA = F^l? in the ab- 
sence of steering, is therefore increased to the larger value 
T — FA^m since the effective capture area of the terrace 
is larger. The additional area for large L is simply dV , 
where the quantity d denotes the in-plane displacement 
of a particle landing on the terrace edge (trajectory t2) 
and V is the island perimeter. This leads to the effective 
flux 



T ^ F{A +dV)^ F{jL^ + eL) 



(16) 



where e — rjd is an adimensional quantity. The actual 
value of d depends on the thermal energy of the incoming 
particles and on the interaction energy between the atom 
and the substrate. For oblique incidence, it also depends 
on the angle of incidence and it should be meant as an 
average over the island perimeter. 

The effect of the parameter e will be illustrated in the 
next Section, in the analysis of experiments on second 
layer nucleation. 

Eq. (|l^) assumes that atoms landing on the terrace 
because they are steered behave as unsteered ones. This 
is not exact, because the extra flux FdV is not deposited 
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FIG. 3: Steering effects for normal deposition. Dotted and 
dashed lines represent the trajectories without and with steer- 
ing, respectively. The adatoms following the trajectories ti 
and t2 reach the terrace only because of the adatom-terrace 
interaction. Adt corresponds to the area inside the dashed- 
double dotted line. 



uniformly, but close to the border of the terrace. Con- 
sequently, the probability W is smaller than the value 
given by Eq. and its use leads to an overestimate of 
steering effects. However, this is a second order effect 
that should be taken into account only if a more rigorous 
theory for deriving the incoming flux J- were employed. 



V. SECOND LAYER NUCLEATION 
EXPERIMENTS 

The nucleation rate lu, as we have discussed in Sec. ||, 
is the rate of dimer formation on top of a terrace of fixed 
size L. Its direct experimental determination is very hard 
because it would require large statistics in single island 
measurements with a fixed size L. 

Let us instead describe how 'second layer nucleation' 
experiments take place. A first possibility (see Sec. V A) 
is to deposit a fraction of a monolayer on the substrate 
and to analyze the statistics of islands with a se cond nu- 
cleus on top. A second possibility (see Sec. VB) consists 
in preparing an ensemble of one monolayer high islands 
which is as uniform as possible, that is with the smallest 
dispersion in size and island-island distance. Afterwards, 
a new fraction of a monolayer is deposited and the statis- 
tics of islands with a second nucleus on top of them is 
studied. 

The relevant quantity is clearly the probability p(t) 
that a nucleation event has occurred on top of an island 
during a time r. For an ensemble of equivalent islands, 
this probability corresponda_to the fraction /(r) of is- 
lands with a second nucleus :E3 



p(t) = fir) = 1 - exp 



dtuj(L{t)) 



(17) 



ond layer nucleation experiment requires two separate 
pieces of information: the nucleation rate lij{L) and the 
growth law L{t) of the islands. The former has a very 
general validity, because it does not depend on the details 
of the experiment. It does not even change if nucleation 
occurs on top of a mound: uj{L) is a 'single-island' prob- 
lem and it is not affected by the growth dynamics of 
the overall surface. The latter piece of information, the 
growth law L{t), is system dependent and an exact cal- 
culation of it is generally impossible: its determination 
involves, in principle, all the surface, in much the same 
way as the problem of submonolayer nucleation does. 

Let us now determine in an approximate way /(t), for 
a uniform array of islands of density A'is [N-^s being the 
number of islands per lattice site). The growth 
is found according to a simple deterministic model£3'E2l 
each island collects the atoms falling in the capture area 
1/Nis- If Ai = jLf is the initial area of the terrace, after 
a time t: 



Ft 

At) = A + j;^^, m 



Ft 



(18) 



Hence, if islands grow from the initial size Li to the final 
size Lf, the fraction of them with a second nucleus on 
top is 



/ = 1 - exp 



dLLu{L) 



(19) 



The integral / dL Luj{L), appearing in Eq. (|T^), can 
not be evaluated analytically for generic ES barriers: the 
reason is the complex L— dependence of W [see Eq. (||)]. 
It is easily evaluated numerically. 

It is useful to define a critical size L^- it is the (final) 
size of the islands, corresponding to a value f — ^. In 
other words, an island grown up to the size Lc has equal 
chance of having or not having a second nucleus on top. 
Lc is determined by the implicit equation: 



dLLuj{L) 



Li 



ln2 F 



(20) 



This expression clearly shows that comparison with a sec- 



We now reanalyze in detail two second layer nucleation 
experiments. 



A. Pt/Pt(lll) 

This experiment is described in Ref. |l3| a nd it has al- 
ready been discussed, according to Eq. (||), in Ref. |l^. 
Using that expression for lo, valid for strong barriers, the 
integral / dL Lto^L) can be easily evaluated analytically. 
Here we do not assume that barriers are strong, we use 
the general expression (^, which will be later shown to 
be equivalent to (Q) , because of the validity of the relation 

< Tdep- 

In the experiment, = and L^, has been evalu- 
ated as follows: after deposition of a dose of platinum. 
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TABLE II: The critical perimeter Vc and the density N^^ of 
islands, for the five different concentrations of CO. 



Pco (10-"mbar) 


0.5 


10 


47 


95 


190 




816 


623 


449 


341 


324 


Nis X 10^ 


2.75 


2.1 


1.83 


1.76 


1.65 



the size L_ of the smahest island with a second nucleus 
on top, and the size of the largest island without a 
second nucleus on top, have been determined experimen- 
tally for several STM topographs. The critical size has 
been approximated bj^he mean of their average values: 

i, = i((i_)+(L+)).y 

Experiments have been performedlla at different CO 
partial pressures, pco- In Table |l| we report Nis and 
Vc for the different values of Pco- 'Pc is the critical 
perimeter and it is three (six) times the critical size Lc, 
for triangular (hexagonal) islands. The values of the 
other physical parameters entering in Eqs. ([5| JT9| , |20| ) are: 
F = 5 X 10~3mL/s, and v = uoexp{-Ed/kBT), with 
1^0 = 5 X IQi^s-i, Ed = 0.2£j2,V, T = 400 K=0.0345 
eV/kg. It is worth remarkingE3E2 that island coalescence 
has already started for the smallest value of pco- So, in 
that case the ES barrier AE is overestimated. 

We now have all the ingredients in order to apply 
Eq. ( po| ) and recover the value of the interlayer hop- 
ping rate, i^', the only unknown quantity. In the hy- 
pothesis of a uniform barrier along the step, v' = 
i'Qexp[— (£'d + /S.E)/ksT\. Two different values will be 
considered for the prefactor, z/q = i/g = 5 x lO^^s""'^ and 
i'q = 2.2 X lO^^s^^. The latter value has been obtained 
by Field Ion Microscopy studies.E3 As a general rule, if 
AE is the ES barrier obtained for a prefactor Vq = vq, 
for a generic value v'q the ES barrier changes to a value 



(A^;)' = AE + k^T\ni,^'Jiyo) 



(21) 



For T ~ 400 K, changing i/'q by a factor two modifies 
the ES barrier by 24 meV; changing Vq by a factor ten 
changes AE by 80 meV. 

In Fig. ^ we report the results for the Ehrlich- 
Schwoebel barrier obtained using the exact nucleation 
rate given in Eq. (P, v'q being equal to the two values 
mentioned more above. An inspection to the original 
STM images (Fig. 1 of Rcf. ^ shows that for the smallest 
and the two largest values of pcoj islands are triangular, 
whilst in the other two cases the shape is less precisely 
defined. Therefore we have considered both hexagonal 
and triangular terraces: the upper and lower values of 
the error bars correspond to such cases. As one can see 
from the figure, at this temperature the effect of the in- 
determination of the terrace shape is comparable to an 
indetermination of Vq of a factor two. 

In Fig. H we compare our results with those obtained 
using the strong barrier approximation, Eq. (^): the lat- 
ter is applicable for AE > 0.2 eV; for smaller values it 
overestimates lu and therefore it underestimates the ES 




le-U le-lO le- 

CO partial pressure (mbar) 

FIG. 4: Values of AE, the ES barrier, for different pco 
values. Circles and diamonds refer, respectively, to i/g = uq = 
5 X lO^^s"-' and to = 2.2 x lO^^s"-'. The error bars are due 
to the indetermination in the island shape. 
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FIG. 5; Results with the exact nucleation rate (for Uq = 
and triangular islands) are compared to Mean Field Theory 
and to the values obtained from the 'strong barrier' approxi- 
mation for oj, see Eq. M). 



barrier. Mean Field Theory, also shown in the same fig- 
ure, is generally incorrect. 

We can now evaluate the ratio between the residence 
time and the deposition time. In the strong barrier limit, 
such ratio is 



= a-f— exp{AE/T)L^ 

Tdep V 



(22) 



For the largest barrier {AE — 0.31 eV), Tres/Tdep < 
0.005. We conclude that for analysing the experiment 
on platinum, the use of formula (^ would be legitimate. 
Notice that, because of the exponential dependence of 
Tres on AE , larger values of the ES barrier quickly inval- 
idate Eq. (0): Tres ^ Tdep for AE = 0.4 eV. 

Let us now consider the effect of steering on the value 
of the ES barrier. We repeat the evaluation of the ES 
barrier by considering an effective flux = F{A + eL) 
for several values of the parameter e, which measures 
the 'strength' of steering. Results are reported in Fig. ^ 
(note the log- linear scale). For triangular terraces on an 
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FIG. 6; The effect of steering on the interpretation of the sec- 
ond layer nucleation experiment on platinum, Ref. Note 
the log-linear scale. Circles and diamonds refer, respectively, 
to pco = 5 X 10~^^mbar and pco = 1-9 x 10~^mbar. Dotted 
lines mark the values of AE in the absence of steering. 



fcc(lll) surface, ry = 3, and e is three times the (maxi- 
mal) in-plane displacement of a particle landing close to 
the terrace edge. It is easy to check that, for the crit- 
ical size of the terrace, the largest value of e (e = 32) 
corresponds to an incoming flux increased by 60% for 
diamonds {pco — 1-9 x lO^^mbar) and by 24% for circles 
{Pco = 5 X lO^^^mbar). The ES barrier is seen to de- 
crease for strong steering, because a larger effective flux 
requires a smaller barrier for producing the same nucle- 
ation rate. The fact that AE is less affected by steering 
when the barrier is smaller is intuitively clear: a smaller 
barrier implies a larger critical size Lc and the larger is 
L, the smaller is the effect of steering, which is an edge- 
effect. 



B. Ag/Ag(lll) 

The experiment on silver reported in Ref. ^ is slightly 
different from the one on platinum. First, a quantity — 
0.1 ML is deposited. Afterwards an equal second dosa 
A9 = 0.1 ML is deposited at two different temperatures^ 
(T = 120, 130 K) and the fraction / of covered islands 
is determined as a function of the island radius r (see 
Table |l|). The density of islands is determined by the 
relation Ni^ — do/Ai = 9o/{'yLf), and AO — 9q implies 
Lf = V2Lj. Furthermore, = 1.1 x lO^^ML/s, i^o = 
2 x 10"s-i and Ed = 0.097 eV. 

Experimental data at different temperatures give the 
possibility — in principle — to determine both the barrier 
AE and the prefactor i/g of the interlayer hopping rate. A 
large debate is going on in the literat| Uijfl a bout the inter- 
pretation of the experimental dataEltScJ and the possi- 
bility that j/Q/i^o 3> 1 (see below). Our contribution is to 
analyze the data with the correct theory and to ascertain 
the possible effect of nonuniform barriers and steering. 

Inserting each pair of experimental values (r, /) in 
Eq. (|l^) yields a value for the ES length, reported in 



TABLE III: Experimental data for the fraction / of covered 
islands as a function of the island radius r, and theoretical 
results for the corresponding ES lengths, £i = lEs{Ti). 





Ti = 120 K 






T2 = 130 K 




r(A) 


/ 




r(A) 


/ 


^2 


7.0 


0.0 




33.5 


0.0345 


5150 


11.5 


0.02 


300 000 


39.0 


0.14 


10 200 


16.5 


0.11 


285 000 


47.0 


0.47 


17000 


21.5 


0.24 


180 000 


55.0 


0.7 


14 800 


27.5 


0.55 


155 000 


63.0 


0.85 


11800 


32.0 


0.88 


200 000 


68.0 


0.93 


11400 


38.0 


1.0 




77.0 


1.0 





Table III, and consequently for the interlayer hopping 
rate, v' ^ through the relation v' — v / i^s- 

For each temperature, we have used the value of f^s 
obtained from the pair (r, /) with / as close as possible 
to / = i, i.e. with r as close as possible to the critical size 
Lc, because it is the statistically most significant point: 
f{Lf = 27.5A)=0.55 at T = Ti, and f{Lf = 47A)=0.47 
at r = T2 . Using the relations 

T2/AT 

, (23) 



AE 



AT \l2 



where AT = T2 - Ti, wc finally obtain AE = 0.30 eV 
and v'q/uq = 1.9 x 10^ 

The value found for the prefactor j/q = 4 x lO^SsT^, 
is huge and it is comparable to a previous analysiaia of 
the same data using the approximate formula (|^) for the 
nucleation rate. 

A possible alternative way to interpret the data is to ar- 
bitrarily set i'q = vq, and allow the ES barrier to depend 
on T. In this way we find AE ~ 0.12 cV at T 120 
K and AE ~ 0.11 eV at T 130 K. Hence we find a 
significant difference for the barrier strength despite two 
rather close temperatures. This difference persists even 
if we suppose that i/q is somewhat larger than vq (e.g., 
i^q/i^o = 10^). The values of AE change according to 
Eq. (pl|), but their difference 

d{AEy = S{AE) - fee AT lii{iy'„/vo) (24) 

is almost unchanged, because k^AT w 10"'^ eV is smaller 
than S{AE) = A£;(120 K)-A£;(130 K)« lO'^ eV. 

We can wonder whether experimental indetermina- 
tions can explain the large value of S{AE), or equiva- 
lently of v'q/vq- As remarked in Ref. ^, island radii have 
been determined experimentally with a precision of about 
SA. Taking into account the indetermination in r, the av- 
erage values of £1,2 are 

« 14 000 ±4 000. (25) 



180 000 ±80 000 



However, even if we consider the smallest £i and the 
largest £2 [in order to minimize Ai?(120 K) and to max- 
imize Ai?(130 K)], we still have that 5{AE) is of order 
0.01 eV [more precisely, S{AE) = 8 meV]. 
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Let us now discuss whether the unexpected experimen- 
tal results can be the effect of non uniform barriers. As- 
suming that the prefactor v'q is equal or even a bit larger 
than i/Q, we have shown that 6{AE) « 10 meV. If a sce- 
nario with two types of barriers applies, Eq. gives 
the relation 



cal integcation is straightforward. In the limit of strong 



barriers! 



L and the analytical integration is possi- 



6(AE) w fceTln ( — 



(26) 



between the difference in the ES barriers and the kink 
concentrations, ci and C2- In the previous equation, T = 
^{Ti+T2) is the average temperature. A value 6{AE) w 
0.01 eV would therefore require that the concentration C2 
of kinks at T = 130 K is two or three times larger than 
the concentratjkin ci of kinks at T = 120 K, which seems 
to be unlikely.Ea 

We finally consider the effect of steering. Increasing the 
value of the parameter e reduces slightly the spreading of 
the ES lengths (third columns in Tab. Ill), but the value 
of i5(A£^) remains practically unaffected. Alternatively, 
if we consider as a free parameter, the ratio Vq/i'o 
is smaller than without steering, but still too large: for 
e = 30 we find AE = 0.23 eV and vy^o = 1.5 x 10^ We 
conclude that neither steering effects nor the presence 
of nonuniform barriers are enough to allow a reasonable 
interpretation of the experimental results of Ref. H. 



VI. DISCUSSION AND CONCLUSIONS 

The present paper has three main goals, i) Provide a 
complete list of the correct formulas to be used for the in- 
terpretation of data from second layer nucleation exper- 
iments; ii) Extend the 'standard' model of irreversible 
nucleation and take into account the effect of steering 
and nonuniform ES barriers; iii) Apply the theoretical 
framework to reconsider some experimental results. Let 
us discuss these issues in detaiLj— 

It isLJiow well establishedQEaEa that Mean Field 
Theorycj is not appropriate to study the problem of nu- 
cleation on top of a terrace jand the reason of its failure 
has been clearly understood^ In Sec. || we give the most 
general formula for the nucleation rate, Eq. ([^), and sev- 
eral approximate formulas which can be used in-.the rel- 
evant limits. With respect to our recent papersa on the 
same problem, we are now able to provide a very good 
simple analytical expression for the probability W that 
two atoms meet (Eq. ^), valid for any barrier strength. 
We also provide the numerical values for all parameters 
appearing in W and in T^es , which depend on the shape 
of the terrace and on the symmetry of the underlying 
lattice, see Table | and App. ^ 

Mainly because of the complicated expression for W , 
it is not possible to integrate analytically the nucleation 
rate, Eqs. ( p^pO| ), and write explicitly, e.g., the criti- 
cal size Lc as a function of all the physical parame- 
ters, F^T,v^v' ^ . . .. This is not a real limitation, in- 
deed, because performing such one dimensional numeri- 



ble. However, that approximation may be inaccurate, as 
we have shown in the case of platinum below AE = 0.2 
eV. 



In Sec. Ill we have considered the possibility that 



the additional ES barrier at steps is not homogeneous. 
The case of a generic distribution of barriers can not be 
tackled analytically, but if barriers are everywhere large 
enough to keep the adatom density uniform, the prob- 
lem is solvable. For a simple bimodal distribution (the 
barrier is equal to AE^ in special kink sites and equal to 
the larger value AEq elsewhere) Eq. (|l^) indicates that 
the system behaves as it had a single effective barrier 
AE^ff equal to the kink barrier plus a correction depend- 
ing on the kink concentration and temperature. This 
result shows that the same experimental system displays 
different effective barriers at different temperatures, or 
at different growth conditions. 

In Sec. 1^ we have studied the effects of the 'incoming 
atom'-substrate interaction (steering) on the total flux 
landing on a top terrace. A detailed treatment would 
require to consider a realistic inter-atomic potential and 
to specify the energy and the inclination of the incoming 
particles. We have introduced a minimal model where 
the 'strength' of steering and all the above variables are 
included in a single adimensional parameter, e = -qd. rj 
is a geometric factor and d represents the (maximal) in- 
plane deviation of an incoming atom. 

Finally, in Sec. we reconsider the data concerning 
two experiments of second layer nucleation, which had 
already been analyzed in Ref. nsL using the strong barrier 
approximation, Eq. (|^). 

In the case of Ft/Ft (111), the additional barrier varies 
from a small value {AE ^0.1 eV) for the clean sample, 
to a value of order 0.3 eV when terrace steps are fully 
decorated with CO. We stress that these values should 
be meant as effective ones. If steering plays a promi- 
nent role, the effective barrier decreases, because more 
adatoms are expected to land on the terrace. If the land- 
ing point of an adatom is displaced up to five interatomic 
distances, d ~ 5, the largest barrier AE — 0.31 eV de- 
creases by 20 meV; ii d — 10, the reduction is twice as 

The case of Ag/Ag(lll) is very debatedElB and the 
question of the actual values of AE and ly' is still open, re- 
quiring additional experiments at different temperatures. 
The interpretation of the data by Bromann et ala can be 
easily summed up. If the prefactor v'q is assumed equal 
to i/Q, the barrier is of order 0.12 eV and it differs by 
about 10 meV at the two temperatures, Ti = 120 K 
and T2 = 130 K. This difference, S{AE), between the 
two barriers can be reduced if i/'q is allowed to increase. 
However, only an exponentially large value of Vq/vq [see 
Eq. (24)] would imply a noticeable reduction of 8{AE). 
If fg/t'o — 2 X 10'', the two barriers are both equal to 
0.30 eV. Data on the fraction of covered islands at dif- 
ferent sizes give different barriers. The spreading of AE 
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at a given temperature can be accounted for by the error 
bar on the determination of the radius r; however, this 
is not enough to explain the difference between AE{Ti) 
and AE{T2). Steering effects are not able to make the 
difference consistent either. In principle, we expect two 
different effective barriers, at the two temperatures. Ap- 
plication of the formula ( p^ suggests that the value of 
6{AE) can be explained by a larger kink concentration 
at the highest temperature, Ck{T2) / Ck(Ti) w 2, but it 
is difficult to find a physical motivation for such a large 
ratio. We conclude pbi^ remarking that a completely dif- 
ferent measurementcH of the interlayer rate in the same 
system, but at a higher temperature, T = 300 K, gives 
a barrier AE = 0.13 eV and a prefactor lApi the same 
order of vq [see also the discussion in Ref . |2J(b)] . 

We believe that possible further developments in the 
problem of nucleation should include the following points: 
a direct quantitative determination of the parameter e; 
an improvement of Eq. to include the case of weak 
barriers; a better assessment of the growth law L{t) for 
the terrace; additional quantitative and controlled sec- 
ond layer nucleation experiments at different tempera- 
tures and for other systems as well. 



the total number of atoms changes by dA = 2jdL L. If 
we set so that only edge sites are comprised in the 
circular ring, we have 77 — 27 d^. Let <f> be the angle 
between the vector joining an edge site with the center 
of the circle, and a nearest neighbor bond. If = 0, 
= 1, whilst for a generic 0, c?l = cos(/). The value of 
dL entering in the relation between 77 and 7 is just the 
average over of cos 0. Because of the different lattice 
symmetries, the average is evaluated between zero and 
7r/4, for a square lattice, and between zero and 7r/6, for 
a triangular lattice. We therefore obtain dL — 2\/2/'k 
in the former case, and d^ — Z/tt 'm the latter case. In 
conclusion, we obtain 77 = 4-\/2 for the (100) lattice and 
ry = 4V3 for the (111) lattice. 

The derivation of Az, the (average) number of missing 
neighbours per edge site, is easily found from V . Let us 
first determine the number 7V„„ of nn bonds which are 
cut by a circle of radius L. It is elementary to write that 
Nnn = 8L for a square lattice and 7V„„ = 8^/3L for a 
triangular lattice. The quantity Az is nothing but the 
ratio between Nnn and V, so that: Az = \/2 for the 
(100) lattice and Az = 2 for the (111) lattice. 



APPENDIX A: INTRA AND INTER-LAYER 
HOPPING RATES 

The quantity v is the total hopping rate on a flat sur- 
face. The jump rate in a given direction is Vd = i^/z^ 
where z is the coordination number (z — A,6 for a square, 
triangular lattice respectively). 

If an atom is on an edge site, its coordination number 
is reduced by Az, the number of missing neighbors. The 
directed hopping rate towards a missing site (i.e. towards 
the lower terrace) is equal to v'j^ and it may differ from 
the hopping rate, Vd, towards a site of the same terrace. 
If we define ly' = zv'^, the absence of additional step-edge 
barriers is equivalent to write v' ~ v. 

Within this formalism, the escape rate from the terrace 
for an atom on an edge site is equal to Az v'^ = {Az/ z)v' . 
Other authors, see e.g. Ref. ||, prefer defining v' as the 
escape rate {v' — Azz/^), but in that case the absence 
of additional ES barriers does not match the condition 
v' = V. 



APPENDIX B: CALCULATION OF SOME 
PARAMETERS FOR A CIRCULAR TERRACE 

The derivation of 7 for a circular terrace is straight- 
forward, because the number A of atoms contained in a 
circle of radius L is equal to the area ttL^ divided by the 
area per atom, A^. Since Aa — \- for a square lattice and 
Aa = •\/3/2 for a triangular one, we have 7 = tt in the 
former case, and 7 — 27r/\/3 in the latter case. 

The derivation of ry, that is, of the number V of edge 
sites, is less trivial. If we change L by a quantity di. 



APPENDIX C: CALCULATION OF A^di. 

We want to evaluate analytically the average number of 
distinct sites N^x^ visited by an adatom during its random 
walk on the terrace. Let us imagine that we perform iV,u„ 
times the following procedure. We let a random walker 
start from a site on the terrace and follow its trajectory 
until it leaves the terrace. For each run r and each terrace 
site s we define the quantity n(r, s) to be 1 if site s has not 
been visited during run r and otherwise. The average 
probability that site s is not visited is 

lun run 

All runs can be grouped according to the number 
F > 1 of 'traversals' across the terrace, defined as fol- 
lows. Once the atom has arrived on an edge site, it 
has the probability poji* ~ Az/z to attempt to move 
outside and pin = 1 — Az/z to move to another site 
of the terrjace: in the former case, the atom has the 
probabilityu a — £es/ (1 + ^es) to stay there and the prob- 
abihty 1 — a = 1/(1 + £es) to leave the terrace. The first 
traversal starts when the atom is deposited on the ter- 
race and terminates when the atom reaches an edge site 
and tries to leave the terrace. At this point the atom 
may leave the terrace, in which case it has performed one 
traversal only, or it may stay on the terrace and start a 
new traversal. According to the this definition, we can 
write 

oc 

n(s)= lim J2N^^^Mr,s))F, (C2) 

run 
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where N^^ is the total number of runs made up of F 
traversals, and < • • • >p is the average value on those 
runs only. 

The quantity n(r, s) is different from zero if and only 
if the site s has not been visited in any traversal. We can 
write explicitly 



n{r, s) = Y\_ nk{rF,s) 



(C3) 



fc=i 



where is the variable n referred to the fc— th traversal 
of the run labelled rp. Inserting ( |C3| ) in Eq. ( |C2| ) we 
obtain 



(C4) 



k=l 



where fi{F) = N^_^^/Nr^^ is the probability that a single 
run is made up of F traversals. 

If a run goes on for F traversals, it means that the 
atom has done F attempts to descend, failing the first 
{F — 1) times and succeeding the last one. Therefore, 



fi{F) = a^-\l - a) 



(C5) 



We now make the approximation that distinct traver- 
sals are independent: 



(Jl nk{rF,s)) ~ {ni{rF,s)) Y[{nk{rF, s)) 



(C6) 



k=l 



k=2 



Using Eqs. (|C5|,|C6|) in (|CJ), we obtain 



n{s) 



F=l 



(1 



'{ni{rF,s))l[{nk{rF,s)) . (C7) 



k=2 



Let us remind the meaning of {nk{rF,s))- It is the 
probability that site s is not visited during the fc— th 
traversal of a run composed of F traversals. This quan- 
tity does not depend on F. We make the approxima- 
tion that it does not depend on fc either, if fc > 2. In 
simple words, we just define (ni(s)) for the first traver- 
sal and {nBs{s)) for all the subsequent traversals. The 



lowerscript BB means that the traversal starts from and 
arrives at the boundary of the terrace. 
Within this approximation, 

oo 

n{s) ~ (ni(s))(l-a)^a^-i(n3B(s))^-i 



F=l 



(1-a) 



1 - a{n-Bj,{s)) 



l+4s(l-(nBB(s))) 



(C8) 



The last approximation is to neglect the s— dependence 
in (nBB(s)) and (rii(s)). Using the simple relation, N^^^ — 
J2s(^ ~ obtain 

A 4s^-+^ ^^^^ 

where N^UL,0) = (1 - {ni{s)))/A and A^r(i) = (1 - 
(?t-bb(s)))/.4 is the number of distinct sites visited by a 
particle starting from the boundary and reaching it again. 
Inserting the expression for N^^{L) = NiL/ \n{L/ Li) we 
obtain the explicit form of W 



w = 



NoL\ti{L/Li) 
HL/Lo) 



Ni£^s+lLlii{L/Li) 



(CIO) 
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